Expression du produit de 2 matrices de Lorentz 
dans le cas general . 


A ° Guimier 


Résume : 
Soit 2 matrices de Lorentz donnée sous leur forme polaire 
M'= A( B').a: etM"= A( B") -Q"', on donne l'expression algébrique 
+ —— 
du produit M = M'° M" en fonction de B', B", Q! Q". 


Un exemple numérique est traite pour illustrer la formule démontrée . 
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Expression du produit de 2 matrices de Lorentz 
dans le cas general : 


Resume : 
_ — | 
Soit 2 matrices de Lorentz M'= Al B').0: et M" = A( B") -Q"'" on donne l'expression algebrique 
> 
du produit M = M'° M" en fonction de B', B", Q', Q". 
Un exemple numérique est traité pour illustrer la formule démontrée. 


On rappelle maintenant qu'une matrice de Lorentz M peut s'écrire comme un produit d 
‘une matrice symétrique À 


et d'une matrice orthogonale @ et cela d'une manière unique : 
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Commençons par étudier le produit de 2 matrices de Lorentz sans rotation : 
Si M=M"'M" est le produit de 2 matrices de Lorentz sans rotation avec : 
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Donc C'!= Id,3 + DBL “14,3 - rl: ceauipermer d'écrire 8 
> > 
en fonction de B', B" : 
> 
Q=|1d4,, - j-PRL. (yyB'B"+c'cr") avec : 
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Considérons maintenant le produit M de 2 matrices de Lorentz quelconques M" et M"' avec 


(16) (y8') 
>) A Q B')o' 
u-=Mm'm"=A(B)Q,m-| ? P , M'= k Lis 
ni — 
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En faisant le produit on voit immédiatement que y= Y'-Y"- ( ‘B'.0':p" +1 ) et que 





tR' + C "Q! [1 , — 
B- = rB’ = p' ce qui donne la partie symétrique de M : A( B) ; 
y (8'a"B" +1) 
Nous allons calculer @ + Pour cela remarquons que : 


0 


Si À est une matrice orthogonale et À = , une matrice orthogonale de Lorentz et 
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t," >) A ‘ > 
(4) A B) A = A( AB) : il suffit de faire le produit matriciel . 


A(B) = 


À 1 0 
On remarque en particulier si on note Q= 0 





A(B).2-0('a A(B).0) -QA('o8) 


Maintenant revenons au produit M : 
— 5 > —, = 
M = M. .M"=A(B').Q AB"). Q''=O'A Ale). A(B"). Q'' 
— — 
or on sait faire le on des 2 matrices symétriques de Lorentz Ale Q'B ) et A( B ") 
— 
en reprenant les formules antérieures, en remplaçant B'par op, 
B par B" et Qpar Q""': 


Al'orB) - A(B") = A(B").or" avec : 





sn. 7 COP) + CP 
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puisque P) =x( 9) D 

On écrit donc Al Q'B' B'). A(B")sous la forme A(°Q Q'B' B'). A( r) = A(B").a"" s 

D'où : 

M= à A(pÉ) Gr dr dr A(pé). (ra). Gr êr- (à. A(pé) à). à.ê.êr 
= A( Q'B ").Q .Q'".0"" et par unicite de la décomposition d'une matrice 

inversible en un produit d'une matrice symétrique et d une matrice orthogonale on a : ni 
M=A(B).Q- A( aB").a0'.0'".0o" onaQ=Q"'.0'".Q"" et par ailleurs que B = Q'B'". 


Donnons maintenant une vérification numérique des formules démontrées : 


#Programme de verification 
#On donne un exemple illustrant la formule demontree . 


#Au préalable on donne des programmes qui nous seront utile par la suite 
restart : 


100 

Id = | 0 1 0 | :unprotect(7y) : 
0 0 1 

f —proc(m) 


LR RS 
return M; 2 m; 3 M; 4 


Ms 2 Myz M, 4 


end proc: 
h :=proc(m) 


1 0 0 0 


0 m;,; M; M,3 














return 
0 m,; M; M; 
0 m;; M; M3 
end proc : 
g :=proc( B) 
local y Bx, By, B : 
1 
7 = = : Bx == BI]: By = BI2]:B = BI3]: 
1—B -B 
} Bx-y Bv-y Bz-Y 
“à ; Fa A 
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Bx-y Hi si PP a 
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end proc: 


#Un proramme qui genére de façon aléatoires des matrices de Lorentz : 
#Cas géneral : 


— 
#On genére aléatoirement le vecteur B ce qui permet de construire la partie symétrique de la matrice . 
— — — 


#On génere aléatoirement trois vecteurs @l, @ et @ et par la methode de Gram — Schmidt 

#on obtient un système ortonormal de 3 vecteurs ce qui permet d'avoir la partie ortnonormale 

#de la matrice. 

#En retour on aura la matrice gerale , vérification qu'elle est bien de Lorentz , la partie symétrique, 
#la partie orthogonale et la vérification que cette dernière est bien orthogonale et finalement B 
#On se limite à 10 chiffres significatifs. 


with(RandomTools) : with(LinearAlgebra) : 
randomize( ) : 


pI =proc() 
local G, co, B 3 M, À Bx, Br, B, ©, GS, œl, @,  : 
unprotect( y) : 
110 0 0 
0 -1 0 0 
G == 
0 0 -I1 0 
0 0 0 -1 


co == evalf( 3 ) ‘ 


Generate( float(range =- 1.1, digits = 10 )) 





co 


Generate( float(range =- 1.1, digits = 10) ) 
co 





Generate( float(range =- 1.1, digits = 10)) 
co 


Bx = BI1]: By = BI2]: Be = BI3]: 





1 
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V1-(B +B +Bé) 


y Bx-y Bv-Y Bz-y 








2 ; 2 
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Bx-y Der or CB Hrr 
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Generate( float(range =- 10.10, digits = 10 
œl :— | Generate( float(range =- 10.10, digits = 10 ; 
Generate( float(range =- 10.10, digits = 10 
float 


( ( )) 
( ( )) 
( ( )) 
Generate( (range =-10..10, digits = 10) ) 
o :— | Generate( float(range =-10..10, digits =10) ) |; 
( ( )) 
( ( )) 
( ( )) 
)) 


Generate( float(range =- 10.10, digits = 10 
Generate( float(range =- 10.10, digits = 10 
@ :— | Generate(float(range =- 10.10, digits = 10 ; 
Generate( float(range =- 10.10, digits = 10 
GS == GramSchmidt( { œ1, @, @ }, normalized) ; 


1 0 0 0 
= 0 GS[1][1] GS[21[1] GS[3][1] 
0 GSI1][2] GS[21[2] GS(3]1[2] | 
0 GS[11[3] GS[2113] GS[3]1[3] 
M = ÀA-Q 
Bx 
return M, MT °G-M, A aao-o, By 
Be 
end proc: 
#Un programme qui donne les parties symétriques et orthogonale d'une matrice de Lorentz 
1 0 0 M{2][2] MI2][3] MI2][4] 
#par la formule :C = || 0 1 0 |- av MI312] M{31131] M[31[4] 
0 01 MI41121 MT4131 MT4114] 


# En entree on a la matricede Lorentz consideree. 


#En retour on a la partie orthogonale Q de la matrice M, la vérification que 


#0 est bien orthogonale et que M = À-Q. 
p4 ==proc(M) 

local B, Bx, By, Bz y Cinv, N, À Q: 

unprotect( 7) : 


y= M1, 1]: 
M£2, 1] 
7 
M13, 1] 


Ciny=|010!-——-p-hB 
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% 
return "M=", M, B, à 0 Q " A: 
end proc: 














#Donnée de M1 et M2 et M = MI - M2 
a = pl():M1 = al1}; fl — af[6]: O1 —al4]: Y —MI[1][1]: 


a == pl() : M2 = a[1]}; B2 — a[6]: 2 —a[4]:2 = M2I1]1[1]: 
M := M1 -M2: 
#Calcul des parametres de M selon 2 manieres : 


7= n-2.(B1" -f(a1) -B2 +1); w=MLU, 1} 


1 7 
C1 = Id + pi - pr”: 
1+Y 
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2 % 
C2 = Id + Ep .pr"": 
l+2 




















M2, 1] 
MIL, 1] 
y - BI + CI -f(01) -B2 M3, 1] 
BE ur BE | y] L 
(BI -f(Qr) -B2+1)-y 
M4, 1] 
ML, 1] 
# Calcul du produit intermédiaire 
%T 
2 %T %T 2 
CI — ne 2 -Utan pr). (tan pr) :C2 = Id + Ep pr", 
1+Y 1+2 
Fe y. (on -pr) + cr. pr 
%T : 
see a 
%T 
anale) one 1 


ICE (on) -pr). B*" + c1 -c2); 





%T 
1+B 

#Calcul de Q 

Q:=f(01) -0 -f(@); 
#Verification de l'orthogonalite de Q : 
0-07: 

#Vérification avec un autre programme : 
Oy = f(p4(M) 141); 


#Calcul de À 
A = g(f(Q1) -B3); 


#Verification avec un autre programme : 


Av := p4(M) [6]; 
#Vérification ultime : 


A := g(f(@1) -B);Mf == A-h(Q);Mv = M; 





1.67468551100000  —0.786693928656870  0.152452215757818  —1.07816628862018 
M 0.847421806900000  —0.262286367684665  0.788973259660187  —1.01333645718720 
7 | 0.543319685500000  0.0222809515380425  —0.627838230770929  —0.948956793240696 
—0.889523221600000  1.24482800374528  0.0811297070426241  0.484843041413859 
1.17560683700000  0.429704069249763  0.368763942601609  —0.247828572708692 
up | —1:109397508100000  0.0281109462715498  —0.795176919184823  —0.615525188553972 
7 | 0.159600840600000  —0.453459541284826  0.557939415794519  —0.713127256086562 
0.587035937400000  0.989055006419444  0.438616143292397  —0.417129788916830 
1.44624323586195  —0.437992209308167  0.855282447038959  0.410212454925417 
0.555984223454041  —1.00324546777367  0.516796223221964  —0.188516984260923 
7 | —0.0209823936898763 —0.419778212709573 —0.583884053574633  0.695202072552982 
—0.884342003856170  0.0955089428148110  —1.05995714088802  —0.805872002215380 
Y = 1.44624323570801 
p = 1.44624323586195 
0.384433413161911 
B—| —0.0145082053021837 
—0.611475291378323 
0.384433413182174 
B —| —0.0145082052379460 
—0.611475291242492 
1.23138695552407  —0.0448401249527581  0.317116487093209 
C1:—| —0.0448401249527581  1.00868949937660  —0.0614535200293941 
0.317116487093209  —0.0614535200293941  1.43460905632545 
1.00550090878963  —0.00802531686643268  —0.0295183245490957 
C2:—| —0.00802531686643268  1.01170819464013  0.0430645038723615 
—0.0295183245490957  0.0430645038723615  1.15839773344865 
B = 1.44624323570801 
By = 1.446243236 
0.975554098856391  —0.0298585577336796 —0.217721534445259 
O3 —| 0.0511654485505307  0.994360729876032  0.0928915267173489 
0.213720136800983  —0.101760529573219  0.971580412367140 
—0.903698229747633  0.322406907778476  —0.281750416398464 
Q=—| —0423535045291322 —0.576547958669593  0.698720626762929 
—0.0628297304816074  —0.750763769356264  —0.657575994735920 





1.00000000177299 
—1.41145609960312 10° 
1.1919624427535210 7° 


0.999999997499976 
2.2447368506190710 710 


—1.4114560996031210 7°  1.1919624427535210 7° 
2.2447374057305910 10 
1.00000000126335 


—0.903698229395731  0.322406908116973  —0.281750416658580 
OQy—| —0.423535044753133  —0.576547958472331  0.698720626679344 
—0.0628297303917543 —0.750763769572174 —0.657575994964109 
1.446243236 0.555984223336941 —0.0209823939197258  —0.884342004369736 
0.555984223336941 1.12636456265454 —0.00476889616615003 —0.200993995744189 
7 —0.0209823939197258 —0.00476889616615003 1.00017997427574 0.00758535047791896 
—0.884342004369736  —0.200993995744189 0.00758535047791896 1.31969869935496 
1.44624323586195 0.555984223454041 —0.0209823936898763  —0.884342003856171 
0.555984223454041 1.12636456268866 —0.00476889611419295 —0.200993995639406 
Éd —0.0209823936898763 —0.00476889611419295 1.00017997427177 0.00758535038927398 
—0.884342003856171 —0.200993995639406 0.00758535038927398 1.31969869893530 
1.446243236 0.555984223336941 —0.0209823939197258  —0.884342004369736 
0.555984223336941 1.12636456265454 —0.00476889616615003 —0.200993995744189 
Fe —0.0209823939197258 —0.00476889616615003 1.00017997427574 0.00758535047791896 
—0.884342004369736  —0.200993995744189 0.00758535047791896 1.31969869935496 
1.44624323600000 —0.437992209449980  0.855282447202797 0.410212455343449 
0.555984223336941 —1.00324546808663  0.516796222889292  —0.188516983975137 
Th —0.0209823939197258 —0.419778213207286 —0.583884053854022  0.695202072596314 
—0.884342004369736  0.0955089427883110  —1.05995714092605  —0.805872002145374 
1.44624323586195 —0.437992209308167  0.855282447038959 0.410212454925417 
0.555984223454041 —1.00324546777367  0.516796223221964 —0.188516984260923 
Lun —0.0209823936898763 —0.419778212709573 —0.583884053574633  0.695202072552982 
—0.884342003856170  0.0955089428148110  —1.05995714088802  —0.805872002215380 
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